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Abstract 

We consider the finite volume approximation of a reaction-diffusion system with fast reversible reaction. 
We deduce from a priori estimates that the approximate solution converges to the weak solution of the 
reaction - diffusion problem and satisfies estimates which do not depend on the chemical kinetics factor. It 
follows that the solution converges to the solution of a nonlinear diffusion problem, as the size of the volume 
elements and the time steps converge to zero while the kinetic rate tends to infinity. 
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1 Introduction 

In this paper, we consider chemical systems with fast reactions where mean reaction times vary from ap- 
proximately 10 -14 second to 1 minute. In particular, reactions that involve bond making or breaking are not 
likely to occur in less than 10~ 13 second. Moreover, chemical systems almost always involve some elementary 
reaction steps that are reversible and fast. 

The study of reactions with rates that are outside of the time frame of ordinary laboratory operations re- 
quires specialized instrumentation, techniques and ways of proceeding (see for example Espenson [1J Chapter 
11]). This work tries to give an efficient, quick and cheap way for numerical investigations of such reactions. 

In this article, we consider a reversible chemical reaction between mobile species A and B, that takes 
place inside a bounded region n C R d where d = 1, 2 or 3. If the region is isolated and diffusion is modelled 
by Fick's law, this leads to the reaction-diffusion system of partial differential equations 

u t = aAu — ak(rA(u) — rs(v)) in Q, x (0,T), 
v t = bAv + f3k(r A (u) ~ r B {v)) in Q x (0,T), 

where T > and Q, is a bounded set of R d . An example of explicit expressions and values for a, (3, k, ta, Tb, a, b 
is given in Section [6] We supplement the system (JTJ by the homogeneous Neumann boundary conditions 

Vu-n = V«-n = on dClx(0,T), (2) 

and the initial conditions of the form 

u(x, 0) = uo{x), v(x,0) = Vq(x) in Q. (3) 

In the sequel we call the system $Q together with the boundary conditions @ and the initial conditions ©, 
Problem T k . 

For a reversible reaction mA f± nB one has a = — m, j3 = n and the rate functions are of the form 
ta(w) = u m and rB(v) = v n . Further discussion about this motivation and some concrete examples can be 
found in Erdi and Toth [9] and Espenson [4]. 

In practice, especially for ionic or radical reactions, changes due to reaction are often very fast compared 
to diffusive effects. This corresponds to a large rate constant k. Bothe and Hilhorst [T] study the limit to 
an instantaneous reaction. They exploit a natural Lyapunov functional and use compactness arguments to 
prove that 

u k ^>u and v k -» v in L 2 (fi x (0, T)) , 
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as k tends to infinity, where (u k ,v k ) is the solution of Problem V k and the limit (u, v) is determined by 



u v 

r A {u) = r B {v) and h - = w, (4) 

a p 



where w is the unique weak solution of the nonlinear diffusion problem 

w t = A</>(w) in !5x(0,T) 
d<f>(w) 



dn 



= on 9Qx(0,T) ^ 



w(x, 0) = wo(x) ■— — uo(x) + ^vo(x) in O, 
a p 



with 



a. , b \ (I. , 1 
— id + t:^ ° — id + — r) , .... 

a 0'J \a P'J (6) 
77 = r^ 1 o r A - 

The identities in © can be explained as follows: the first one states that the system is in chemical equilibrium, 
while the second one defines w as the quantity that is conserved under the chemical reaction. Given a 
function w, the system ((4} can be uniquely solved for (u, v) if ta, tb are strictly increasing with for instance 
tvl(R + ) C r_g(R + ) so that r\ = r^ 1 or a is well defined and strictly increasing. Under these assumptions u is 
the unique solution of 

1 1 , s 

~ u + -T, r i\ u ) = w, 

a p 

which gives the explicit representation of u and v 

1., 1 \-\ ^ 1 



We assume the following hypotheses, which we denote by TL: 

1. Let Q, be an open, connected and bounded subset of M. d , where d — 1,2 or 3, with a smooth boundary 

on, 

2. uo(x), vq(x) G L°°(Q) and there exist constants U, V > such that < uo(x) < U and < Vo(x) < V 
in fi, 

3. a, f3, a, b and k are strictly positive real values (sometimes we use the notation ka — a and kj3 = j3), 

4- Let rA(x), rs(x) G C 1 (R) be strictly increasing functions, such that rA{0) = rs(0) = 0, and assume 
furthermore that r A (R + ) C r B (E+). 

We recall from Bothe and Hilhorst [T] Section 2] that Problem V k has a unique classical solution (u k ,v k ) 
on every finite time interval [0, T], for all nonnegative bounded initial data. By classical solution, we mean a 
function pair (u k , v k ) such that u k , v k G C 2A (n x (0,T]j n C 1 ' (Ti x (0,T]) with u k ,v k G C([0,T]; L 2 (fi)) 
(see also Ladyzenskaja, Solonnikov and Ural'ceva [llj). 

Next we present a notion of a weak solution of Problem V k , which will be used in the sections [4] and \5\ 

Definition 1.1. We say that (u k ,v k ) is a weak solution to Problem V k if and only if 

1. u k ,v k G L 2 (0,T;H 1 (n)) and u k ,v k G L 2 (0,T; 

2. Let \& be the set of test functions, defined as 

* = [V 1 £ C 2,1 (n x [0,T]) : W n = on 9nx[0,T] and %j}{T) = o}. 
For a.e. t G (0, T) and all Tp £ $ 

/ tto(x)i/'(a;, 0) d:r + / u k (x,t)^t(x,t) dx + a / u k (x, t)Atj)(x, t) dx 
Jn Jn Jn 

-aj ip{x,t)[r A {u k {x,t)) ~r B (v k (x,t))"j dx = (8) 
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and 



I Vo(x)ip(x, 0) dx + / v k (x,t)ip t (x,t) dx + b I u k (x,t)Aip(x,t) dx 
Jn Jn Jn 

+ $ [ iP(x,t)(r A {u k (x,t)) -r B (v k (x,t))^ dx = 0. (9) 



We remark that every essentially bounded weak solution of Problem V k , in the sense of Definition 11.11 
is also a classical solution. 

This paper is organized as follows. In section 2 we define a finite volume discretisation and an approxi- 
mate solution for Problem V k ■ In section 3 we prove a discrete comparison principle which yields 
discrete L°° estimates, and we show the existence and uniqueness of the approximate solution. Section 
4 contains technical lemmas used further in the convergence proofs. The convergence of the approximate 
solution to the classical solution of Problem V k in the case of fixed k is proved in section 5. In section 6 we 
use a suitable Lyapunov function and we obtain a discrete L 2 (0, T; H 1 (Q)^ estimate, which does not depend 
on k. We then apply Kolmogorov's theorem and deduce the convergence of the approximate solutions to 
the classical solution of Problem V k . Afterwards we show that the approximate solution {u^jVj)) converges 
to (u, v) defined in (J7]| as k tends to oo and the size of the discretisation parameters tends to zero. 
In Section [6] we present numerical results obtained with our finite volume scheme, for the reversible dimeri- 
sation of o-phenylenedioxydimethylsilane (2, 2-dimethyl-l, 2, 3-benzodioxasilole) which is a reaction of the 
type 2A B (see Meyer, Klein and Weiss 12 ). On the one hand we compute the approximate solution 
{ut},Vx>) of the solution (u ,v ) of Problem V k and on the other hand the numerical approximation w-d of 
the solution w of the problem and we check that 

for k large enough and size (T>) small enough. 

Remark 1.2. In what follows we denote by C , Ck and positive generic constants which may vary from 
line to line. 



2 The finite volume scheme 

The finite volume method has first been developed by engineers in order to study complex, coupled physical 
problems where the conservation of quantities such as masses, energy or impulsion must be carefully respected 
by the approximate solution. Another advantage of this method is that a large variety of meshes can be used 
in the computations. The finite volume methods are particularly well suited for numerical investigations of 
conservations laws. They are one of the most popular methods among the engineers performing computations 
for industrial purposes: the modelling of flows in porous media, problems related to oil recovery, questions 
related to hydrology, such as the numerical approximation of a stationary incompressible Navier - Stokes 
equations. 

For a comprehensive discussion about the finite volume method, we refer to Eymard, Gallouet and Herbin 
[B] and the references therein. 

Following [B], we define a finite volume discretization of Qt- 

Definition 2.1 (Admissible mesh of Q). An admissible mesh M of Q. is given by a set of open, bounded 
subsets of Q. (control volumes) and a family of points (one per control volume), satisfying the following 
properties 

1. The closure of the union of all the control volumes is Q.. We denote by tuk the measure of each volume 
element K and 

size (M) — max rrih- 

K£M 

2. K n L = for any (K, L) G M 2 , such that K ^ L. If K n L ^ 0, then it is a subset of a hyperplane 
in R d . Let us denote by £ C T 2 the set of pairs (K, L) , such that the d — 1 Lebesgue measure of K n L 
is strictly positive. For (K, L) £ £ we write K\L for the set K n L and m K \ L for the d — 1 Lebesgue 
measure of K\L. 

3. For any K £ M we also define Mk = {L £ T, (K,L) G £} and assume that dK = K\K = 

(Znan)u(u^^li). 
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4- There exists a family of points (xk)k£M, such that xk £ K and if L £ Mk then the straight line 
(xk,%l) is orthogonal to K\L. We set 

di<\L = d(x K ,x L ) and T K \ L = A|L , 

<j>k\l 

where the last quantity is sometimes called the transmissibility across the edge K\L. 

Since Problem Problem V is a time evolution problem, we also need to discretize the time interval (0,T). 

Definition 2.2 (Time discretization). A time discretization of the interval (0,T) is given by an integer 
value N and by a strictly increasing sequence of real values (t' n ^) n g{o,...,JV+i} with t^ ' = and v- N+1 > = T, 
The time steps are defined by 

t (n) = t (n+l) _ t (n) fgr „ £ {Q, . . . , N}. 



We may then define a discretization of the whole domain Qt in the following way. 
Definition 2.3 (Discretization of Qt)- A finite volume discretization T> of Qt is defined as 



e{o,...,JV+i} I ! 



where M, £ and (xk)kgt are given in Definition ] 2.1] and the sequence (t < -"-') n g{o,...,JV+i} * s a time dis- 
cretization of (0,T) in the sense of Definition ] 2. 2] One then sets 



i (V) = max {size (M),^ : n £{0,...,iV}}. 



We present below the finite volume scheme which we use and define approximate solutions. We assume 
that the hypotheses ft hold and suppose that T> be an admissible discretization of Qt in the sense of 
Definition 12.31 We prescribe the approximate initial conditions 

u)9 = — — / uo(x) dx and v^P = — — / vo(x)dx. 
m K J K m K J K 

where K £ M, and associate to Problem V k the finite volume scheme 

m K {v% +1) - uP) - tf\ Tkh(u^ +1) ~ <4" +1) ) + t™m K a(r A (u% +1) ) 

m K (v^ +1) - vP) - t™b £ T K[L (v£ +1) ~ - t^m K p{r A {u^ +1) ) 

L&M K 

Note that is a nonlinear system of equations in the unknowns 

/ (n + l) (n+l)\ 

\ U K > V K JKEM, ne{0,...,JV}- 

For and t £ (0, T) let K € M be such that x £ K and n £ {0, . . . , N} be such that t m = 0, t^ N+1) = T 

and t £ (tS n ^ , t' n+1 '] . We can then define the approximate solutions 

ut>(x, t) — u^ +1 ' > and Vr>(x,t) = v^ +1 ^ . (12) 

In the next section, we will prove the existence and uniqueness of the solution of the discrete problem (llip . 
together with the initial values (|10|l . 



(10) 
(11) 

-r B (^ +1) ))=0. 



3 The approximate solution 

In this section we prove the existence and uniqueness of the solution of the system . Let us start with a 
discrete version of the comparison principle. 
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Proposition 3.1 (Discrete comparison principle). We suppose that the hypotheses TL are satisfied. Let T> 
be a discretization as in Definition \2.3[ Let and (u^\v^)kgm be given sequences of real 

values such that 



u$ ^ and Vjp ^ 



for all K S M. If the sequences {u^ +1 \v^ +1 ^)k^m, ne{o,...,iV} an d 

\ v£ +V> )kg.M, n6{o,...,JV} satisfy the equations (fTTj) with the initial values (u^ ,v^)ksm and (u^ ,v^)} 
respectively, then for K £ M and n G {0, . . . , TV} 

u^'af 1 ' and v% +1) <v% +1) . (13) 

Proof We set = uP - u^ l) and v^ = - for all K G M and n G {0, . . . , TV + 1} and define 



^ 1) =(rx(«J +1) )-rA(fiir 1, ))/^ 1) : 



whenever u K 7^ (else A K = 0) or v K 7^ (else B K ' = 0). Since the functions rA and tb are 
monotone increasing, it follows that A£ and 73^ +1 ' are nonnegative. We then have, by subtracting the 
discrete equation (|11[) for u£ +V> and for u^ +1 \ 



i^if 



+ t} ' m K a[ r B {v K K j - r B {v y K 
for K £ M. and n G {0, . . . , TV}. Setting s + = max(s, 0) and using that s ^ s + , (s + t) + ^ s + +t + we obtain 



(15) 



^mK(4 n) ) + + 4 n) « E ^i44 n+1) ) + +4 n) ^"(^K n+1) )-^(4"- 

where K £ M and n G {0, . . . , TV}. Next we multiply the inequality (|15|l by indicator of the set where 

is nonnegative. Since the right-hand-side of (|15p is nonnegative as well, we obtain, acting similarly for both 

components, 

M'l.. r" 1 _ t. . \ \(u (n+i ^ 



mJl + t^UA^ + ^- E 



, /~(n)\ + . An) rri i~(n + l)\ + , An) ~( I (n + l)\ f ~{n+Vj 

^ Wk(%') +*j « 2^ T k\l{uI ') +t\ 'mKa[r B {v} ( ')-r B (v K K 
(«)/' 3 £>(«+!) , ± T 7- \ \ ,„<■" I '-M 



(16) 



<mK(4 n) ) + +4' l) k E T , if|L(fii' l+1) ) + +4 n) ^A'/3(r A (n< f " +1) )-^(4" 



,-,("+!) 



L£Af K 



Since 



b% +i) (* ( k +i y =(^(4 n+i) ) -r B («r i) )) + 

we add the first equation of (|16p divided by a and the second equation of (|16[) divided by f3, which yields 



<w,ifT-«) ^ 2^ t a'|lK J + +^ -5 2^ T K\L(vi ') , (17) 
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for K £ M. and n 6 {0, . . . , N}. Let us note that 



E E 



KeM L£Af K 



- e e T L]K (^ +i) r 

E E t k]l (^ +1) ) + 



LeM K<=M L 

Summing the inequalities (I17|l over A" € M, we get 



E E ^(4" +1) y 

KM £€A/" K 



which therefore leads, by induction, to 



E 



inn 



.(n+l)s + 



1 



'K 



0. 



where n £ {0, . . . , N}. It implies that + = {v^^ L ') + = 0, which completes the proof. 

Corollary 3.2 (Discrete contraction in L 1 property). With the notation from Provosition \3. li we have that 



E mK 



(n + l) -(n + l)| 
ML — Mi, 



+ 



| (n + l) -(n + l) 







E 



m K 



(n) -(n)| 
\ U k ~ U k \ 



+ 



I («) 

Vi - VI 



p 



for n £ {0, . . . ,N}. In other words, the discrete counterpart of the L 1 (Q.) -contraction property for solutions 
of (fTJ) (see e.g. f7^) is preserved by the numerical scheme (| 1 1[) . 



Proof The proof directly follows from the proof of Proposition 13.11 Let us consider the term uk- We 



multiply the equation (I14|l by sgn 



Then, the inequality x ^ \x\ yields 



m K 



I -(n+l) I 



«S m K \u y K '\ +t\ 2^ T k\l\u\ >\ +4 >m K a\r B (v y K >)~r B {v y K ') 



We proceed in the same way for v 



(n+l) 



and remark that 



A 



B 



(n+l) I -(n+l) I 



= \r A (u^>)-r A {u^)\ 



(n+l) I, (n+l) | 



J x ) 



= \rB(v K )-r B {' 



.(n+l) 



which enable us to obtain the counterpart of the inequalities in (|17p which we sum over K £ M, as in the 
proof of Proposition 13.11 This yields the result. ■ 

We now are in a position to prove a discrete L°° estimate for the approximate solution. 

Theorem 3.3. Let T> = (^J^A,V,£, (4"' ) )ne{o,...,]v+i}) be an admissible discretization of Qt in the sense of 

i be given by (| 10|) and 



DefiniUon \2.3l We suppose that the hypotheses Ti are satisfied. Let (u^ 



(u 



(n+l) ,,("+1) 
K ' V K 



satisfy (JTTJ) for K £ M and n £ {0, . . . , N}. Then 

u ( " +1) ^U+^V and SC v ( " +1) ^ V + @-U, 



(18) 



for all K £ A4 and n £ {0, . . . , iV}, where U and V are the positive constants from the hypothesis Ti.2. 

Proof From Proposition 13 . 1 1 we immediately obtain that u^ +V) and Vj£ are nonnegative for K £ M and 
n £ {0, . . . , iV}. In order to find a discrete upper solution, we consider approximate solutions of the corre- 
sponding system of ordinary differential equations. More precisely, we consider sequences (M' n ')ng{o,...,iv+i}> 
(i)' n ') n g{o jv+i} (we postpone for a moment the proof that they exist) such that 



[/, v 



(<>) 



V 



G 



and 

Q (»+i) _ G C») = afct (») ( rB ( € («+D) _ r^(u ( " +1) )) , 



-fn + 1) 

IT ' — IT ' 



/3fc4" ) (r A (x t (" +1 >)-r s («(" +1 ))), 



(19) 



for n £ {0, . . . ,7V}. We note that the sequences (u' n+1 - > ) ng { jy}, («' n+1 ^)ne{o,...,JV} satisfy (|11[) with the 

initial data U,V. Therefore they satisfy the comparison principle from Proposition [XT] which yields 

sC u {n+1) and ^ « (n+1) for all n 6 {0, . . . , TV}. (20) 

Adding up the first equation of (|19[) divided by a and the second one divided by (3, we obtain 

u^+V _U V 

a (3 a (3 a (3 

We deduce from the previous equation and from (|20p that 

0<u (n+1) ^U+^-V and 0sJ« (n+1) ^V+@-U, (21) 

p a 

for n £ {0, ...,7V}. ■ 

In order to prove the existence of the sequences (u'"') ng { ,...,]v+i} an d (i' l - n ') ne {o,...,JV+i} we use the 
topological degree theory in finite dimensional spaces. The reader can find basic definitions as well as 
further informations about this powerful theory in Deimling [3]. An example of the application of this tool 
to the analysis of finite volume schemes can be found in Eymard, Gallouet, Ghilani and Herbin [5]. 

With T, g : R 2 -> R 2 defined as 

0(u< n >,€ (B) ) = foi'"-"^!'"')-^!.'"))), 
we rewrite the system (|19[) in the form 

T/x-(n+l) -fn + l)\ | n / -ln + \) -(n+X)\ f-(n) -<n)\ 

J-(u K , ir ' ) + yhr , tr ') = y := [u y , tr 
Moreover we see that setting O = 73(0, r) C R 2 a ball centered at (0,0) with a radius 



r > 



yj(u + %v) 2 + (v + £uy 



we fulfill all the assumptions of Theorem [21 Theorem 3.1, page 16]. 
For the continuous function Tt : [0, 1] x 73(0, r) — > R 2 given by 

u/\ -(n+i) r-(n+ll\ t-i „-(n.+l) -(n + l)\ i \sif-(n+i) -(n + l)\ 

we see that 

d(T+ AS, 73, (u (n) ,v {n) )) = d(H(X),B, (u (n) , «<">)) 

for all A £ [0, 1] and (u (n) ), (v (n) ) such that < M (n) < U - fV and < t) (n) < V - g (7 for n € {0, . . . , TV}. 
On the other hand we deduce from [3] Theorem 3.1, page 16] (dl) that 

d(H(0),B,{u {n) ,v (n) )) = 1. (22) 

In view of [3] Theorem 3.1, page 16] (d3) and (d4), (122[) implies that the equality 

Fivr ',v K ) + <5(tr ',V S ) — (it ',IT ') 



has a solution or, in other words, that there exists a solution of JTSJ. The uniqueness of this solution 
immediately follows from Proposition 13.11 
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We can prove in the same way the existence and uniqueness of the solution of the system (|lip . Indeed, 
we rewrite in the form 

t({u^ +1) ) k ^{v^ +1) ) k&m ) +G((u% +1) ) KeM ,(v% +1) ) KeM ) = ((uP) KeM ,(vP) KeM ), (23) 

where T, Q : R 2e — > R 2e , with the number of control volumes for the discretization V, are continuous 
functions given by 

G{{u^) K& M,{^K)KeM) = {m,m), 

where 

K LeU K 

and where 

We set O = B(0, R) C R 2e a ball centered at zero with a radius 

i ?> ^/e((/ + |\/) 2 + e(i/ + |(/) 2 . 

Since > 1, we deduce from the discrete L cc (Qt) estimate of Theorem [33] that the equation (|23[l does not 
have any solutions on dO. Applying again !3, Theorem 3.1, page 16] with TL{\) = T + XQ and A 6 [0, 1] 
completes the proof of the following result. 

Theorem 3.4. We suppose that the hypotheses TL are satisfied. Let T> be a discretization as in Definition 
\2.'A Let (u^p ,v^P)k£M be given by (JTUJ). Then there exists one and only one sequence 

, (n + l) (>i + l)\ 

(U K ,V K )K£M, n€{0 JV}i 

which satisfies (| 1 1 p . with the initial condition (u^ , v^)k eM ■ H 



4 Convergence proof with k fixed 

We begin with the discrete version of L 2 (Qt) estimates of the gradient of the approximate solutions. 

Proposition 4.1. We suppose that the hypotheses Ti are satisfied. Let T> be a discretization as in Definition 
\2.3l Let ((TO]) and (fTTjl give the sequences (uk ,v^')keM and ( u k +1 \ v k +1) )k£M, ne{o,...,N}, respectively. 
Then, there exists a constant Ct > 0, which does not depend on T>, but which depend on all the data of the 
continuous Problem V k (namely, the constants a, (3, U, V including k and the functions va, tb), such that 



N 

(n+l) „(n+l)> 



KeM KeM n-0 {K,L)e£ 



(24) 



and 



5 E ™«K +1) ) a -5 E ™*(<4 0) ) 2 +*E4 n) E t^v^-v^Y^cu 

(25) 



KeM KeM n=o (K,L)ee 



Proof For the sake of simplicity we only present the proof for the M-component. We multiply the first 
equation in the finite volume scheme (|11[) by u^ +1 ^ and sum the result over all K £ M and over all 
n £ {0, . . . , N} to obtain 

Si + S 2 + S s = 0, (26) 
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where 



<si = E E rnK { u k +1 



) „.(»)K>+ 1 ) 



n=0 KeM 



* 2 = -«E4 n) E E M M 

n=0 K6.M L€Af K 



K U K > 



(n+1) + (n + 1) 



i V 

53 = «E4 n) E 



n=0 -K"6A1 



Since 



we deduce that 



(4 i+i) ) 2 - = \{u^ +i) ) 2 - ) 2+ K u - +1) - 



* = iE E — ((«i? +1, ) a -(«$?) j + ^E E 

N / \ 

^EE-(K n+1) ) 2 "K n) ) 2 )- ^ 



n=0 KeA<l 



All the terms in the sum on n on the right hand side of (|27p simplify except for the first and the last ones. 
We have that 

Si > \ E — (4f +1) ) 2 - 5 E m K (u^y. (28) 



KG.M JfGA-f 



We can perform a discrete integration by parts to obtain 

N 

<s 2 = «E4 n) E T K|i (4" +1) -4 n+1) ) ■ (29) 

n=0 (K,L)e£ 



Finally we use the hypothesis Ti 4 and the inequalities in (118[l to estimate the last term, namely 

N _ 

-S, E ^4U+- V)+r B {V + ^U))(U + ^V) 

N 

s5afcCE4 n) E ™x = afcC|«|T, (30) 

n=0 A'6A4 

with some positive constant C. 

Identities (|26[) and (|29[) together with the inequalities (128p and pop immediately give (|24[) . Since the 
argument in the case of the ^-component is similar, we omit the proof. ■ 



4.1 Space and time translates of approximate solutions 

We now turn to the space translates estimates. We use here methods which have been presented for example 
by Eymard, Gutnic and Hilhorst [8] and by Eymard, Gallouet, Hilhorst and Slimane [7J. The results of 
the current and the next subsection together with the technical Proposition 17.21 will imply the relative 
compactness of the sequence of approximate solutions. 

Proposition 4.2 (Space translates estimates). We assume that 

1. T> = ^A4 , £, (x k) k gT , (^"■ ) ) )1 e{o,...,jv+i}) is an admissible discretization of Qt in the sense of Defini- 
tion \2.S\ 

2. the hypotheses Ti and assumptions 1)361) are satisfied, 

3. the functions [ut> and v-p) are derived from the scheme (|10p - (|lip and given by the formulas (11211 . 
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Then there exists a positive constant Ck, which does not depend on D, but depends on all the data of the 
continuous Problem V k , including k, such that 

f f (uv(x + £,t) - UT ,(x,t)) 2 dxdt 4 C k \£\(2 size {!>) + |£|), (31) 
o Ja^ 

and 

[ [ (vv(x + £,t)-v-p(x,t)) 2 dxdt^ C fc |£|(2size(D)) + |£|), (32) 

Jo J(l£ 

for all £ € R d and for Q% defined as in Proposition \7.S\ 

Proof Inequalities (|31[1 and (|32|) follow from the estimates (|24p and (|25|l. respectively. We refer to [6j 
Lemma 3.3] for a details. ■ 

Proposition 4.3 (Time translates estimates). Let the assumptions of Proposition \4~1^ are satisfied. Then, 
there exists some constant Ck > 0, which does not depend on T>, but which depend on all the data including 
k, such that 

/ (ur>(x,t + t) — UT>{x,t)) 2 dx dt < Cfc(size (T>) + r) (33) 

J(JX(0,T-t) 



and 



Jn-, 



(vT>(x,t + T)-v v {x,i)) dx dt < C k (size (V) + r) , (34) 

/f2x(0,T-r) 

for allr e (0,T). 

Proof In order to apply Lemma |7. II (see appendix), we follow the same steps as in [8] Lemma 5.5]. The 
only difference appear in the nonlinear part of the equations. However, these can be easily estimated using 
the regularity properties of functions r^(-) and rs(-), as well as -L°° estimates (|18|l in Theorem[ 



4.2 Convergence proof 

In this section, we state convergence results with k fixed. This differs from next section where we will 
introduce additional hypotheses about the nonlinear reaction terms and obtain convergence results which 
permit us to pass to the limit as k — > oo. 

Theorem 4.4. We suppose that the hypotheses TC are satisfied. Let (ut>,vt>) be the approximate solution 
defined by (1101) . (|lll) and (|12[1 . There exist a pair of functions (u k ,v k ) and a sequence (ux> m , wz> m )meN of 
(ut>,vt>) such that 

(wp m ,«z> m )meN converges to (u k , v k ) € (L 2 (0, T; ff 1 ^))) 2 

strongly in L 2 (Qt) as size(X> m ) tends to zero. The function pair (u k ,v k ) is a weak solution of Problem V h 
in the sense of Definition \l.l\ 

Since Problem V k is a uniformly parabolic system, (u k ,v k ) must coincide with the unique classical solution 
of Problem V k ■ This immediately yields the following result. 

Corollary 4.5. The pair (ut>,vt>) converges to the unique classical solution (u k ,v k ) of Problem V k as 
size (T>) tends to zero. ■ 

Proof of Theorem 14.41 In view of the estimates (|3ip . (|33[l and Proposition 17. 21 which is a consequence of 
the Frechet-Kolmogorov Theorem [2J Theorem IV. 25, page 72], we deduce the relative compactness of the 
set («d) so that there exists a sequence of (ui> m )^_ 1 and a function Uk, such that uu m — ► Uk strongly in 
L 2 (Q T ) and weakly in L 2 (0,T; H 1 (Q)) , as m -> oo. 

The same conclusion holds for the v-component. Indeed, the inequalities (|32j) and (I34p permit to apply 
the compactness result in Proposition 17.21 for the sequence {v-D m ) rn _ 1 - There exists a function Vfc such that 
Vu m — » Vk strongly in L 2 (Qt) and weakly in L 2 (0, T; _H" 1 (J1)) as m — » oo. 

Next we show that (Uk, Vk) is a weak solution of Problem V , in the sense of Definition ll.il Since the 
proof for the ^-component is similar, we only present here the detailed proof in case of the it-component. 
Let ip £ ty, where ^ is the class of test functions from Definition ll.il We multiply the first equation of (|lip 
by tp(xK,t^ n '), where tp € Then we sum over all K £ M and n G {0 . . . N — 1} to obtain 

7-U /-T-U I t~l — II / i 

Ira — I 2m + 2 3m — u , 
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where 



T 1-=E E ^(4 n+1) -4"V(«:* W ). 
JV-l 

r 2 " m = a £ t<»> ]T £ r K|£ (ui™ +1 ' - w £ +1) M**,* (n) )> 



n=0 KeMLeU K 
JV-l 



7L = E4"' E ^«(^(4 n+1) )-^(^ n+1) ))v'(^,t (n) ). 



n=0 KS.M 

The complete proof, that 



/ uo(x)4>(x,0) dx — / / Uk(x,t)ip t (x,t) dxdt 
Jn Jo Jn 



and 



lim — ~ a / / Uk{x,t)Aip(x,t)dxdt, 
m -*°° Jo Jn 

can by found in [8] Lemma 5.5]. Let us focus on the proof that 



lim T£ m = ak / (rA(£4) — r B (Vk))i>(x,t) dx dt. 



m — >oo 



We write 



JV-l 

& 

r-T 

a 



E4' l) E m^rx^J-rB^ 1 '))^^,*"") 
ip(x,t) (r A (Uk) - rs(Vfc)) dx dt 

/o 

EE/ / (rA(u^)-r B (v^))^(x K ,t^)dxdt 

n=0 KeM JtM Jk 
JV-l /■i<' l + 1 ) /. 

-«E E / / ^(^MWO-r^Vfc^dzdt 

-a V / / ip(x,t)(r A (Uk) ~r B (V k ))dxdt. 

Thanks to the regularity of the function ip, the last sum above converges to zero . Moreover, 

jv-i r t < - n+1 '> r 

"EE/ / (rA(u% +1) )-r B (v£ +1) ))^(x K ,t^)dxdt 

n =OK£M Jt(n) Jk 

JV-l - t (»+l) „ 

""EE/ / ^(M)MZ4)-r s (V fe ))dxdi 

JV-l . t («+l) „ 

=«E E / / (v-(^,t (n) )-^,t))x 



n=0 ' 



x (r A (u^ +1) )~r B (v^ +1) )dxdt 

JV-l „*(" + !) „ 

+ «E E / / ^(M)Mu£ +1) )-r-A(Z4))dzd* 

„_n a,, JtM J K 



JV-l / , t (" + 1 ) r 

"EE/ / ^(^i)(r S (^ +l3 )-r S (V fc ))d^dt 



(35) 



Next we show that the three terms above tend to zero asm-t oo. First we take their absolute value and 
apply the triangle inequality. The Cauchy-Schwarz inequality applied to the first sum of the right hand side 
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of yields 



EE / / k(xK,t^)-^(x,t)\\rA(u% +1) )-r B (v% +1) \dxdt 

N-l . ( ("+U r v 1/2 

EE/ / (lK*K,t (n) )-^M)) a da!dt) x 

n=0 Jt( " ) Jk / 

/W-l / .i<" + 1 > /■ 

HEE/ / M4 n+1 VM^ +1) f 



"i E U / ' f'- JM - •*^ :l) ') 2 d: ' : 

n = if £.M ' 

The first term of above product converges to zero, as m — * oo, since ip(x,t) is smooth enough. The second 
term is bounded. Indeed, it is sufficient to remark that ta{UjI ), Tb{Vk ) are bounded for all K £ M 
and n £ {0, . . . , iV}. The last two terms in (|35p are similar and we show how the proof goes with the first 
one. Indeed, 

"EE/ / Wx,t)\\rA(u% +1) ) - r A (U k )\dxdt 



n=0 if 6A-1 * ( K 

JV-i / .tC"+ 1 ' 



|rA(w^ +1) )-r A (W fc )|dxdt 



< a|MU-<o T ) EE(, / 

n=0 ifeAl^* 1 "' ' /A ' 
iV-1 „ t (" + l) - 

^ a||^lk~(Q T ) E E / / ll r Alli~(QT)l'"( x K.* ( " ) ) - dxdt 

n.OJf£M JlW 7K 

/JV-l ,t< n+1 > „ vl/2 

< a||^IU«»(Qr)( EE/ / ll r Alli°°(Q T ) d:rdt ) x 

V„ = 0KM Jtl " ) Jk / 

JV-l /.*(" + !) f. xl/2 

EE/ / h(^A-,t (n) )-w fc (x,t)| 2 d^dt . 



n=0 if GA4 ' 

The L°° norm of the function ta(x) is taken over the finite interval [0, U + %V\. Because ta{x) is of class 
C 1 (72.) the first term of the above product is bounded. The second one converges to zero since u-p m —* Uk 
as m — * oo in L 2 (Qt)- B 

5 The case that k tends to infinity 

In order to prove the convergence of the finite volume scheme when size (73) tends to zero and k tends to 
infinity, we impose some additional conditions on the nonlinear terms va(x) and tb{x). At first we prove a 
counterpart of Proposition 14.11 

Proposition 5.1. Let us assume hypotheses 7i. Moreover, we assume that the functions rA(x), tb(x) 
satisfy 

r K e C X (R), r' K {-) > on (0,+oo), 

rm n a v sr '^ ^ (36) 
r K (U) — 0, and nmsup < oo, 

where k £ {A, B}. Let T) = LM, £, [xk)k^t, (t^)«e{o jv+i}) ^e an admissible discretization of Qt 

in the sense of Definition \2.3i and the sequences , VjP)iceM an d ( u i? +1 '; w ic +1 ') ksm ne{o n} are 
given by (|I0l) and (|If p . respectively. Then, there exists some positive constant C which is independent of 
the discretization T> and of the reaction rate k, such that 

E4"'( E T m (u^-u^y+ £ T K]L (v^-^ +1) f)<C (37) 

1=0 ^(K,L)€E {K,L)e£ ' 

and 

N 

fc E4 n) E m K (r A (u^)~r B (v^)) 2 < C. (38) 

n=0 K£M 
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Remark 5.2. Observe that the condition (I36|) holds, for example, in the case that the rate functions ta{s) 
and 7\b(s) behave like s 1 , for some positive 7, whenever s — » + . 

Proof of Proposition I5TT1 Let (a,b) 6 (R + ) 2 be such that rA(a) = rs(6). We define two functions 

r fl (6) 7 S r s (a) 

which are continuous on R + because of hypotheses H and the assumptions (|36|1 . We can extend these 
functions to also be continuous at s = 0. To do so for the function Va{s) we check that the hypotheses (|36[) 
sr a ( s] 

give the integrability of £-r- on the interval [0, a] and we pass to the limit 

rA(s) 

,. T7 , , h 1 f o-r' A {a) , 1 sr' A (s) 

lim Vk(s) = - / — da hm s < oo ; 

s^o+ a J rA(cr) a s _o+ rA(s) 

where we have applied de l'Hospital theorem as formulated in [131 Theorem 2, p. 174]. 
For a given e £ (0, 1) and n 6 {0, . . . , iV} we consider 



Since Vl'(s) = > for all s £ R+ we deduce that 

ar A (s) 

V A (u^ +e) - Va{u^ 1] + s) = ( U W - 4" +1) )^(4 n+1) + e) + ^ - 4" +1) )M'(*) 

> (^-4 n+1) )Vl(^ +1) +s) ■ 

As a consequence 

We substitute mg- — itjf ) from the scheme (|lip . which yields 

sC 4 n+1) (e) + 4 n+1) (e), 



A (n + x) (£) = _ t wa ^ +e ) -yi(4 n+1) + e )} 



with 

(if,L)e£ 
and 



Since there exists some constant C > such that Va'(s) ^ C for all s G [0, £/ + §V] we can use the L°° 
bound (|18[) and the mean value theorem to obtain 

4" +1) ( £ K-4 n) ^ E r K|L K" +1) -4 l+1) ) 2 - 

(K,£)e£ 

Following the same steps for the function 

fl C™+D (£)= £ mx(^(4" +1) +e)-^(4 re) +^), 

we arrive at 

A (n+1) (e) + B (n+X) (e) ^ _ £,(»+!) _ D (n+1) ^ ( g g) 

with 



(K,L)ee (k,l)£S 



„(«+i) _„,(™+ih 2 
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and 

J D ( " +1) (e) = t^k J2 m K (r A (u^ +1) ) - r B (v£ +1) )) ( lnr A (^ +1) + e) - \nr B {v { £ +1) + e)), 
where we used that aV A (s) = In and that (3V B (s) = In - rs ' 



Let 

D[ n+1) (e) = i< n) fc E m K (r A {u { Z +1) + e) - r B (t4" +1) + e)) (lnr A (^ +1) +e) - lnr s (t4" +1) +e)). 

ifSAI 

Since the inequality 

(c-d)(lnc-W) > i£_^L 
c + a 

holds for all (c, d) £ (o, U + ^ V + l) x (o, V + ~U + l) , then 

D[ n+1) (e) > tWi J]) m ff ( r ,( u r+^^(r , +^)) 2 - 

6 if £M 

where Ct is an upper bound for r A {c) + r B (d) with (c,d) G (0,17+^V + l) X (0,V+-U + l). Such bounds 

p a 

exist in view of Theorem 13.31 and regularity of the functions r A (x) and r B (x). Let us define 
B (n+i) (e) :=4«)A53 mK ( rA(u ^) +E) _ rfl(t7 ^) +e) )». 



The assumptions (|36|) and Lemma 17.31 in the appendix, imply that for all K £ M and e > small enough, 
there exist constants Ci, C2 such that 

Mwk" +1) +e) -r A {u% +1) )\ < Cie, |r s (4" +1) + e) - r B {v% +1) )\ < C l£j 

and 

|lnr A (u^ +1) +e)| < C 2 (| lne| + l), | lnr s («^ +1) + e)| < C 2 (| lne| + l). 

Then 

_ D (n+i) ^ Ce \n\ k (\ ln£ | + _ D («+D( e ) ; 
for some positive constant C. As a consequence 

A (n+1) (e) + B {n+1) (e) < - C (n+1) - D (n+1) (s) 

< - C (n+1) - £ (n+1) (e) + Ge|n|fc(| lne| + l). 
Now it is possible to pass to the limit in (|39p . We obtain 

A (n+1) (0) + S (n+1) (0) < -C (n+1) -£ (n+1) (0). 

which is 

E rn K (v4 U ^)-V A (u^))+ E m K (v B (v^) - V B (vP)) 

+4 n) ^| E ^„K" +1) -4r i) ) 2 +4 n) ^ E W4 n+1) -<4" +1) ) 2 

(K,L)ee (K,L)es 

E ^(r A (^ +1) ) -r S (4" +1) )) 2 < 0. 
Now we sum the above inequalities over n 6 {0, . . . , N} to obtain 



n=0 (K,i)6£ 



+ 1 ) „,(™+ 1 )\ 2 



>, W 1 N 

n=0 (K,L)ee b n=0 KeAl 



KE.M 
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Since Va(s) and Vb(s) are nonnegative and continuous on [0, +00), and since u^ +1 \ v^ +1 ^ are nonnegative 
and bounded for all K £ M and n € {0, . . . , N}, this concludes the proof. ■ 

5.1 Space and time translates of the approximate solutions 

Since we have already presented the general methods in section [4] we only give here some essential ideas, 
leaving out the details of the proofs. 

We begin with a counterpart of Proposition 14.21 

Proposition 5.3 (Space translates estimates). Let us assume that 

1. ^A^ , £, (xk)k£T, (£ ( '"' > )ng{o,...,jv+i}) is an admissible discretization of Qt in the sense of Definition 

2. hypotheses H and assumptions (|36p are satisfied, 

3. functions (ut> and vt>) are derived from the scheme (110[) - and given by the formulas (|12[1 . 
Then there exists a positive constant C which is independent of T> and k, such that 

[ [ (uT>(x + £,t) -UT>(x,t)) 2 dxdt < C|C|(2size(©) + (40) 

Jo io ( 

and 

[ [ {vT,{x + £„t)-VT,{x,t)) 2 dx dt^C\£,\{2 817,0 (V) + (41) 

Jo Jn^ 

for all £ £ R d and for fij defined as in Proposition \7. 2\ 

Proof As it was in the proof of Proposition 14.21 we refer to [6] [Lemma 3.3] for a complete proof. The only 
difference is to apply the result of Proposition 15.11 instead of that of Proposition 14.11 ■ 



Let us now prove an analogue of Proposition 14.31 
Lemma 5.4 (Time translates estimate). Let the assumptions 1. 2. and 3. of Lemma \5.3\ be satisfied. Set 

wt> = —u-d + -up. Then there exists a positive constant C which is independent ofV and k, such that 

a p 

(wT>(x,t + r) — wi)(x,t)) 2 dx dt ^ C(size (T>) + r), 

(1X(0,T-t) 

for allr E (0,T). 

Proof The proof is similar to that of Proposition 14. 3 1 We present below the essential steps of the argument. 
We define 

A(t) := / (wx>(x,t + t) — wx>(x,t)) 2 dx, 
Jn 

which can be easily transformed into 

A(t)= E ^R" (i+T)+1) -^ (,)+1) ) a 

KeM 

. n(t + r) 

= E (K B(t+T,+1) -«& (t)+1) ) E mx(«^-«w; 

keM ^ n=n(t)+l 



Since 

1 ( (n+l 



(n+l) (n) _ J- / (n+l-) „»\ , 1 ((n + l) (n) 
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we can apply discrete integration by parts in the scheme (|lip to obtain 



n(t + r) 

A(t) = ^ 4 n) E rA 1 LK" +1) -uk" +1) )K n(t+T,+1) -4 n(t+T)+1) ) 

n=ra(t)+l (if,£)S£ 
n(t+r) 

E 4" 3 E rA-,44 n+1) -4 n+1) )K n(t)+1) -^ (t)+1) ) 

n=n(t)+l (K,L)SS 



X? 4" 5 E r^(4 n+1) -<4 n+1) )(* 



- V* t {n) V* T < n+1 ) _ „(™+ 1 )V„„("( t + T )+ 1 ) _ .„(™( t + T )+ 1 ) > \ 

/3 

n=n(t)+l (K,L)g£ 
n(t+r) 

+| E 4" 3 E r K| ^«r i )-«r i) )(4 n(t)+1) --r t,+1) ). 

n=n(t)+l (K,L)eS 

Next we estimate the second term in the sum above, to obtain 

n(i + r) 



~ E 4 n) E \^K n+1) -4 n+1) )-x^iK nW+1) -»r )+1) ) 

n=n(t) + l (K.L)eS 

n(t + r) 

£ E *f° E ^(4 

n=n(t) + l (K,L)eS 



n(t + r) 

a ^ (n) T / „„(n(t)+l) (n(t)+l)^2 



2a 

n=n(t) + l (X,L)G£ 



n(i + r) 

2q 



E 4 B) E T KlL (u^-^ +1) r 



n=n(t) + l (K,L)e£ 
n(t + r) 



a v-^ t , „(™(*)+i) , \ 2 



n=n(t) + l (K,L)e£ 
n(t+r) 

a ^ *00 T , f„Wt)+i) „,(«W+i)\ 2 

rc=n(t)+l (it,L)e£ 

where the first inequality follows from the relation ±siS2 ^ ^( s i + s 2) an d the second one follows from the 

simple inequality (si + S2) 2 ^ 2(si + si)- To conclude the proof we integrate above inequalities over K with 
respect to the time variable t. Next we apply Proposition 17.11 ffor details see [5] Lemma 5.5]). ■ 

Proposition 15.31 together with Proposition 15.41 immediately give the following corollary. 

Corollary 5.5. Let the assumptions 1. 2. and 3. of Proposition^^ be satisfied. We set w-p = —u-d + —vt>- 

a p 

Then, there exists a constant C > 0, which is independent of the discretization parameters T> and of k, such 
that 

f (wv(x + £,t) -WT>(x,t)) 2 dx dt < C|C|(2size(©) + |£|), (42) 

Jn ( x(o,t) 

for all £ G R d and f2j = {x G R d , [x, x + f] C fi}. Moreover 

/ (wT>(x,t + t) - WT>(x,t))' 2 dxdt < C(size(£>) +r), (43) 

J(lX(0,T-r) 

wfeere r G (0,T). ■ 



5.2 The limit as size{V) tends to zero and tends to infinity 

We state below the main convergence results of this paper, first only letting the size of the volume elements 
and the time steps tend to zero, and then also letting the kinetic rate tend to infinity. 



1G 



Proposition 5.6. We suppose that the hypotheses Ti are satisfied. Let k > be arbitrary and let (u k , v k ) £ 
(^L 2 (0,T; be the unique classical solution of Problem V k . Then the sequence (it|>, of approx- 

imate solutions of Problem V k given by (1 lOfl . and (j 12|) converges strongly tn (l 2 (Qt)\ to (u k , v k ) as 



L 



size CD) — > 0. In particular, the sequence — —u^, H v^, converges strongly in L 2 (Qt) to the function 

a p 

w k = — u k + — v k as size CD) —* and there exist positive constants Ci, Ci which do not depend on k, such 
a p 

that 

{r A (u k )-r B (v k )) 2 dx < Ci, (44) 

«x(0,T) 

(w k (x + £,t) ~w k {x,t)f dxdt < C* 2 |£| 2 , (45) 

Q 5 x(0,T) 

(w k (x,t + r) - w k (x,t)) 2 dx dt < C 2 t, (46) 

/f2x(0, T-t) 

where r 6 (0,T), (6l ti and fi 5 = {x G R d , [sr, x + £] C fi}. 

Proof To prove the result we use Corollary 15,51 The method of proof is similar to that of Theorem 14.41 
■ 

It is now possible to pass to the limit as k — > +oo. 

Theorem 5.7. Let (u|>, Vx>) be the sequence of approximate solutions of Problem V k , defined by (|10[1 . 
and GU). Then 



and 



Mb — > ( ^ ld + JjV ) ('"') 



as size CD) — > and k tends to oo, where r\ — r^orA and where w is the unique weak solution of the problem 
©-©. 

Proof Let Wx> = — uv + t; v t>- The estimates from Corollary 15.51 which are uniform with respect to 
a p 

k, permit to apply Proposition 17.21 As a consequence we deduce the relative compactness in L 2 (Qt) of 
the sequence {ujp}. Then there exist a function w £ L 2 (Qt) and a subsequence {w!^ } such that 
converges to w strongly in L 2 (Qt) as ki tends to infinity and sizeCD m ) tends to zero. Theorem [3]3] implies 



that w is nonnegative and bounded in Qt- The inequality ((38}, namely 

ki\\rA(u^ m ) - r B (v^ m )\\ 2 L2{QT) < C 

where the positive constant C is independent of hi and sizeCDm), implies that 

r A-( u -D m ) ~ r s("o m ) ~* in L 2 (Q T ), 
and consequently almost everywhere, as ki tends to infinity. Then 

v v m = r l\ u v m ) + e^ m , 

where n(s) = i'b 1 i rA ( s )) and tends to zero almost everywhere as sizeCDm) tends to zero and ki tends 
to infinity. In view of the hypotheses 7i4 the function r/(s) is well defined on [0, oo). Moreover, 

H(v^ ) = tUp — -je^ — > w a.e. in Qt, 

P 

where H(s) = — s + \r](s). Hypotheses 7i4 ensures that the function H(s) has the continuous inverse 
a p 

function. Then Lebesgue's dominated convergence theorem implies that there exists a function pair (u, v) £ 
(L 2 (Qt)) 2 

Ujj — > u and Djj — > w in i (Qt) (47) 
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as size(T>m) tends to zero and ki tends to infinity. 

Next we identify the limit pair (u, v). Let w be the solution of the problem (J5j ™ ([6]l and the functions u and 
v are denned as in (0, namely 

Let e > be arbitrary. We have that 

ll W ~~ U Wl2(Q t ) ^ \\ U ~ U *\\l*(Q t ) + II" * - U T) m \\l^(Q t ) + \\ U -V m - U \\ L Z{Q T ) ■ 

From 1|47[) we deduce that there exists ko and So such that, for all ki > ko and all size(Dm) < So, the last 
term of the inequality above is less then e/3. From Theorem 1 in [T] there exists some ko, such that for all 
ki > ko, \\u — u ki \\ L 2(q ) ^ e /3- Then, fixing ki = max(fco, ko), we can take size(V m ) < So small enough so 

that by Proposition 15.61 ||n fc ' — it^j | L 2 ( q ) ^ e /3- 

Since the argument for the w-component is similar, this completes the proof. ■ 



6 Numerical example 

In this section we give an example of an application of the finite volume scheme f) in one space dimension. 
For the numerical experiments we choose the reaction of the reversible dimerisation of o-phenylenedioxydi- 
methylsilane (2, 2-dimethyl-l, 2, 3-benzodioxasilole) which has been studied by H NMR spectroscopy. The 
kinetics of this reaction can be described quantitatively by a bimolecular lO-ring formation reaction and a 
monomolecular backreaction (for further details we refer to Meyer, Klein and Weiss |12j). 
Since the reaction is of the type 2A =± B, the reaction terms take the form 

rA{u) = fciii 2 and tb(v) = k^v. 

Moreover a = 2 and = 1. For this particular process benzene was chosen as a solvent. Then it was possible 
to estimate rate constants for both reactions at the temperature T — 298K, 



ki w 1, 072 • I0" 4 L 2 mol" 2 and k 2 ~ 2, 363 ■ 10" 6 L 2 mor 



and diffusion coefficients 



awl, 579 • 10" 9 m 2 s _1 and 6 » 1, 042 • 10" W 



In the experiment we set k = 1 for the chemical kinetics factor. We remark that it is equivalent to the 
situation when coefficients a, b, ki and k^ are of order 1 and k is of order 10 4 . In fact, we can multiply the 
system |T]) by 10 9 and change the time scale as t i— > W 9 t. However the above reasoning is formally correct 
and shows in an explicit way the order of the kinetics factor k; in our example we decided to keep constants 
in the form given by the spectroscopic analysis. 

Figure 1 shows the initial conditions uo(x) and vo(x), defined as 

'0 for x £ [0, 0.03] 

U0 ^ =, i * sin (^(z -0.03)) for x £ [0.03, 0.1], (4§) 



and 



i C os \mx 



wo. 7 for x£ [0,0.07], 

v {x) = < V / (49) 

[0 for x € [0.07, 0.1] 

On Figures 3 and 5 we see the time evolution of the solution (u%,,v^) until the times T^ax = 10 5 [s] and 
Tmax = 10 11 [s], respectively. Then we follow the evolution of the solution w-p of the nonlinear diffusion 
problem ©-([l]) for initial condition deduced from that used in the reaction - diffusion Problem V k ■ We 
have used a uniform mesh with h — 0, 002 and initial time step t$ = 10 -8 to obtain the approximate solution 
(uj>,i>!,) and ts = 10 -6 to obtain the approximate solution wnto the nonlinear diffusion problem. 

We can use the approximate solution wt> and the formulas (0 to define functions iij and v%. Indeed, 

let 

Mb = h(wv) and v% = g(h(wv)), 
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where 



and 



uf \ 1 ( \ 2 , 4fc 2 ak 2 



Q /JK2 



Proceeding in the similar way as in the proof of Theorem 15,71 we write for the u-component that 
||u£> - K^)\\ L 2 {Qt) < || tA - u\\ l2{Qt) 

+ h k - /i W|| L 2 (Qt) + H&H - K™-d)\\ l2{QtY 



(50) 
(51) 

(52) 



We simultaneously pass to the limit as size (T>) — > and k — > oo. From Theorem 1 5. 7 1 we immediately deduce 
that the first term on the right hand side of ()52p tends to zero. The same conclusion holds for the two other 
terms. Indeed, [T| Theorem 1, Sec. 3] implies that j|it fe — /i(tu)|| i2 \ tends to zero as k — > oo. Moreover, 

Theorem 5.1] yields that || w — wt>\ \ ,2 in ~* as s i ze C^) — * an d since the function h is well defined 
and continuous, we conclude that for every small e > there exist T> small enough and k large enough so 
that 

\\ui,-h(wp)\\ L2{QT) <e. 

We proceed in the same way to show that for every small e > there exist T> small enough and k large 
enough so that 

||*4 ~ 9(H u >t>))\\ L 2 (Qt) < £• 
The results from our numerical experiment agree with above analysis, since 



whereas 



max \ (%, T"m; 



- h(wT>)(x,T^)\ 



max Wc^Tmax) - g{h( w v)( x , T mwd)\ 



max Uj>(x, 



ax J 



:4,74- 10 , 
~ 4,032 ■ 10" 

• 3,121 ■ 10" 14 



max !)b(i,t1 x ) - g(h(w-p){x,T^ &K ))\ ~ 1,84- 10" 

E£Si I loo 





0.00 0.01 0.02 0.03 0.04 0.05 0.O6 0.07 



0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 



Figure 1: Initial conditions no and vo, defined in Figure 2: Initial condition wq — jj-Uo + jj v o, where 
and (O. uo, v are defined in ((35J) and 
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Figure 3: Evolution of M ap and v ap on the time Figure 4: Evolution of w ap on the time interval 
interval (0, . . . , T^J. (0, . . . , T^ ax ). 






0.00 0.01 0.02 0.03 0.04 0.05 0.06 



0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 



Figure 5: Evolution of it ap and v ap on the time 
interval (0, . . ..T^J. 



Figure 6: Evolution of ui ap on the time interval 
(0, . . . , T max ) . 



In order to show the accuracy of our method in the fast reaction limit, let us increase the kinetics 
parameter k in Problem V h keeping all other data as previously. Let 



Ju = max \uT>{x,T max ) - h(w-v){x, T max ) 



and 



j£ = max \vi(x, Tmax) 



ff(/i(u>x>)(a:,T I L x )) 



Table 1 below shows the numerical results. 



k 


J u 




O v 




HT 7 


2,3498-10" 


2 


4,5272.10" 


2 


1CT 6 


3,976-10" 


3 


1,988.10" 


3 


10 -5 


4,699-10" 


8 


2,3498.10" 


8 


1CT 4 


9,3725-10" 


10 


4,6862.10" 


10 


icn 3 


2,498-10" 


16 


9,992.10" 


16 


icn 2 


1,7892-10" 


12 


8,95.10" 


13 


io- 1 


4, 163-10" 


17 


9,992.10" 


16 


10° 


3,121-10" 


17 


1,84.10" 


16 



Table 1: The accuracy of our method in the fast reaction limit, when the kinetics parameter k in Problem 

V h increases and all other data are unchanged. 
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7 Appendix 



The proof of the following result can be found in [8j Lemma 5.3 and Lemma 5.4]. 

Lemma 7.1. Let (t'"') n ez be a strictly increasing sequence of real numbers such that lim n ^_oo f'™' = — oo 
and lim n _ 00 t^ = oo. Moreover, let tg := t^ n+1 ^ — t^ be uniformly bounded. For allt£~M> we denote by 
n(t) an integer n, such that t 6 [t^ n \t^ n+1 '). Let (a'™')«gz be a family of nonnegative real values such that 
a («) _^ o j or finitely many n G Z. Then, for all r G (0, +oo) and (el 

/ E (4 n) « (n+1) )* = -E(4 n) « (n+1) ). (53) 

, n(t + r) . 

/( E 4"> (n(t+C)+1) *^(- + max4"))E(4 n) « ( " +1) )- (54) 

The following proposition is a direct corollary from Frechet-Kolmogorov Theorem (see Brezis [3] Theorem 
IV.25, page 72]). 

Proposition 7.2. Let O be a bounded and open subset o/R , d = 1, 2 or 3. £e< (u! n )n.eN ^ e ffl sequence 
of functions w„(x,t) : R d x R — > R, smc/i £/iai 

1. /or aH n G N, w n S L°°(0) and i/iere exists a constant C'b > which does not depend on n, such that 

H^n || ^ Cbj 

2. i/iere eaxsi positive constants Ci, C2 and a sequence of nonnegative real values (jln.)neN, such that 
lim„^oo /j, n — and 

(lU n (x + f,t + r) - W n (x,t)) 2 dx dt < Cl|f I + C2T + 

/or £ G R d , r G R, n G {0, . . . , JV} and 

0^ T = {(x,t) €R d+1 : the interval [(x, t), (x + £, t + t)] hesinO). 

Then there exists a subsequence of (w„)„gn, denoted again by (w n )neH and a function w G L°°(0) such that 
w n — > w in L 2 (0), as n — > 00. 

Lemma 7.3. Le£ A > and a function r(s) satisfying 

r G C 1 (R), r'(-) > on (0,+oo), 

sr'(s) (55) 
r(0) = 0, and limsup — -^—^ < 00 
s ^o+ r ( s ) 

be given. 

Then there exists £0 > on/?/ depending on r(s) and C > 0, depending only on r(s) and A, such that for all 
£ G (0, £0) and /or aiZ w G [0, A] the inequality 



lnr(« + e)| <; C(|lne| + l) 



holds. 



sr f ( si 

Proof Let a = limsup — r-V- + 1. There exists a constant £0 > 0, such that for all s G (0,£o) 
s ^o+ r(s) 

^ < (56) 
r(s) s 

Let e G (0, £0) and u G [0, A]. Then, either r(u + e) ^ 1, which implies 

lnr(u + £) < lnr(A + e ) 

or r(u + e) ^ 1. In that case 

I lnr(u + e) < | lnr(e)|, 
holds. Integrating inequality (|56[) over the interval [£,£0], we obtain 

I lnr(e)| < | lnr(e )| +a(|lne | + |lne|). 

Setting C = max {a, a\ In £0] + | In r(eo)|, In (r(A + £0)) } we conclude the proof. ■ 
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